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Abstract 

Many molecular theories of nematic liquid crystals consider the constituent 
molecules as cylindrically symmetric. In many cases, this approximation may 
be useful. However the molecules of real nematics have lower symmetry. 
Therefore a theory was developed (Mol. Phys. 30 (1975) 1345) for an en- 
semble of such particles based upon a general expansion of the pairwise in- 
termolecular potential together with the molecular field approximation. In 
this study, we would like to handle this molecular field theory by using Tsallis 
thermostatistics which has been commonly used for a decade to study the 
physical systems. With this aim, we would like to investigate the dependence 
of the order parameters on temperature and would like to report the variation 
of the critical values of the order parameters at the transition temperature 
with the entropic index. 
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I. INTRODUCTION 



It is well known that an essential characteristic of compounds forming liquid crystals 
is the rod-like shape of their constituent molecules, with an high length to breadth ratio. 
Therefore the molecules are supposed to be cylindrically symmetrical. For example, the or- 
dering matrix which is often used to describe the partial alignment in a mesophase, contains 
only one independent element and this can be determined by some techniques [1] . The fact 
that the molecular cylindrical symmetry is assumed is appealing to a statistical mechani- 
cian, because the pairwise anisotropic intermolecular potential required in any calculation 
is simple for such particles [2]. However the molecules, in fact, are lath-like and thus do not 
possess the high symmetry. 

The ordering matrix has two principal components and therefore these components are 
required to describe the orientational order of a uniaxial mesophase composed of lath-like 
molecules. In this sense, the deviation of this ordering matrix from cylindrical symmetry 
was found to be significant [3] . The importance of deviations from cylindrical symmetry may 
be inferred from unambiguous determinations of the ordering matrix for rod-like molecules, 
such as 4, 4' — dichlorobiphenyl [4]. Moreover it is found that these matrices are comparable 
to those estimated for a pure mesophase [3] . 

There are some studies in which the consequences of deviations from molecular cyhndrical 
symmetry are investigated. It is shown that a system consisting of particles with a lower 
symmetry than D^oh is capable of existing either as a uniaxial or a biaxial liquid crystal 
[5]. The possible existence of a biaxial phase is studied in detail for a system of hard 
rectangular plates using a lattice model [6], the Landau approach [7] and the molecular 
field approximation [8] . The deviations of the ordering matrix from cylindrical symmetry is 
clearly determined by the molecular symmetry and the element of the ordering matrix for 
the long axis will also be influenced by the form of the pseudo-intermolecular potential. The 
calculations of the ordering matrix for an ensemble of hard rectangular particles is performed 
in [9]. It must be emphasized that although these calculations are of some interest, they may 
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not be particularly realistic because some experiments indicate that dispersion forces may 
make a dominant contribution to the anisotropic intermolecular potential [9]. Considering 
the cases above, Luckhurst et al. developed a theory [10] for non-cylindrically symmetric 
molecules interacting via a completely general intermolecular potential within molecular 
field approximation. 

For a decade, nonextensive statistics has an increasing interest and recently Tsallis ther- 
mostatistics (TT) has been applied to the nematic-isotropic transition [11-13] as a nonex- 
tensive statistics. In [11], the Maier-Saupe mean field theory has been generahzed within 
TT and applied to a nematic liquid crystal, para-azoxyanisole. In the other study, [12], the 
the effects of the nonextensivity on the dimerization process has been studied, and finally 
the mean field theory of anisotropic potentail of rank L — A has been generalized within 
TT and the effect of the nonextensivity on the order parameters has been illustrated in 
[13]. Up to now, the mean field theories for uniaxial nematogens formed by cylindrically 
symmetric molecules have been studied by using TT. In this manner, we aim, in this study, 
to enlarge the applications of TT to the liquid crystal systems and to handle Luckhurst 
et al.'s theory which considers the molecules to be non-cylindrically symmetric. In doing 
so, we first give some essential properties of Luckhurst et al.'s theory. Then we mention 
on TT and its axioms. Finally, we apply TT to the Luckhurst et al.'s theory and some 
possible concluding remarks are made. We must emphasize that we would like to give only 
the possible contributions of the nonextensivity to the theory. So we must keep in mind that 
since one relies on the generalized theory or not, more extensional studies related with it 
must be performed in the nematic-isotropic transition. However, we believe that this study 
is sufficient to give motivation for further apphcations of TT to the liquid crystals. 

A. The Molecular Theory for Uniaxial Nematics Formed by Non-cylindrically 

Symmetric Molecules 

The intermolecular potential for particles of general shape is given by [10] 
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in a product basis of Wigner rotation matrix [14] , where ri2 is the distance between molecules 
1 and 2. The orientation of molecule Qi-r in a coordinate system containing the intermolec- 
ular vector as the z axis is denoted by fli-R. This potential energy is invariant under rotation 
of the coordinate system about z axis. Therefore the summation in Eq.(l) can be restricted 
as follows [10]: 

Uu = J]«L,L.^nin.(ri2) ^L^^U^l-^) ^^^^ • (2) 

In what follows, the redundant subscripts on the coefficient u{ri2) will be suppressed. Be- 
cause it is convenient to define the molecular orientation in terms of a common coordinate 
system, the potential energy Uu could be transformed to a laboratory frame. The choice of 
this coordinate system is determined according to the symmetry of the liquid crystal phase, 
so for a uniaxial mesophase, the laboratory z axis can be taken to be parallel to the sym- 
metry axis of the mesophase. The transformation of U12 is carried out by performing the 
rotation from the intermolecular vector to the molecule in two steps using the relationship 

= E D^^i^i-R), (3) 
j 

where the subscript R — L is the rotation from the laboratory to the intermolecular frame, 
1 — L denotes that from the laboratory to the molecule coordinate system. Then the inter- 
molecular potential can be written as 

C/12 = ^Si(^^i-l) £'S,(^2-l) D^^^Xi^R-L) Dl^^rjnn_L). (4) 

If the distribution function for the intermolecular vector is independent of orientation, then 

one could use the orthogonality of the rotation matrices to evaluate the ensemble average: 

^n,m {^^R-l) Vj^ _^{^Ir_l) 2L1 + I ■ ^ 

The partially averaged potential may then be written as 
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Then now it is the time to average over the orientations adopted by particle 2. However 



one needs D^:^l^^. This average is taken to be independent of the orientation of molecule 1 
within the molecular field approximation and it is to be identified with the normal ensemble 



average. Since we only consider a uniaxial mcsophase, -D^''„2 vanishes when j 7^ and L is 
even. In this manner, the potential is reduced to the following form: 

/IN"! 

L,even 

where the summation is over L ^ (L even). The scalar contribution could be ignored 
because we would only like to study the orientational properties of the mesophase. The 
average of the orientational pseudo-potential for molecule 1 is given by 

^1 = E |^^^-n,n./^&\(fil-L)<l. (8) 
L,even 

The expansion coefficients in Eq.(8) are defined by 

_ j ^LnunnAi'n) <^{-VoiJl■■■r^)/VT) dri...dr/v 
''^™- /exp(-[/„(ri...r^)/A;r) dri...dr^ ' 

where Uo{ji...Y]^) is the scalar potential of the ensemble [10,15]. 
The pseudo-potential is conveniently written as 



um^ E ^^.-pOSjI/^t), (10) 

L,even 



where 



CL.P=E(-iri^' (11) 

m 

with (3 and 7 is the Euler angles and define the orientation of the director in the molecular 
coordinate system. 

Luckhurst et al. made some comments in [10] on this form of the pseudo-potential and 
the symmetry properties of the coefficients CLjp- Because Luckhurst et al. investigated the 



influence of deviations from molecular cylindrical symmetry on various orientational order 
parameters for a uniaxial nematic mesophase, they tried to minimize the number of the 
variables in the calculation without any loss of the essential physics. As a first approximation, 
they considered only those terms with L equal to 2 in the expansion of the pseudo-potential. 
The neglect of terms higher than quadratic could be a good approximation, since similar 
assumptions for cylindrically symmetric molecules had let to results in reasonable accord 
with experiment [10,16]. As mentioned in [10], the number of expansion coefficients are 
reduced by appealing to some specific model for the interactions or by imposing symmetry 
restrictions on the molecules. If one follows Straley [9] and takes the molecules to be hard 
rectangular parallelopipeds, then the coefficients C2jp are independent of the sign of either p 
or j and zero if either of these subscripts is odd, that is 

C200 = {W^ + L^) - 2W (L^ + B^)+L {W^ + B^) + 8WBL\ /3, (12) 

C220 = {L^ - BW) {B - W) (13) 
C222 ^-L{W-Bf/2, (14) 

where L is the length, B the breadth and W the width of the molecule. It must be noted 
that this parametrization is only approximate [9] and also anisotropic repulsive forces are 
probably not dominant in determining the behaviour of real nematics. Therefore assuming 
the molecules to interact via dispersion forces, the expressions for C2jp could be derived. 

The similar restriction can be imposed on the expansion coefficients C2jp, without ap- 
pealing to specific forms of the intermolecular interactions, by using more formal arguments 
based on the molecular symmetry and its influence on the pair potential. For example, if 
each molecule has a centre of symmetry, then [10] 

(ri2) 

= KL,L„-n2(n2) (15) 

and consequently 

CLjp = (-1)^' c,-jp = (-1)^ c,j-p = (-1)^'+^ c^-,-p. (16) 
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In addition, if the molecules also possess a plane of symmetry orthogonal to their z axes, 
then 



ULiL2mnin2{l^l2) — ( — 1)'^^ 1tLiL2m-nin2 (^12) — 1)'^^ ^L^Li 



mni —712 



[1-12) 



^LiL2m— m— n2 



(17) 



and 



CLjp = (-1)^ CL-jp = (-1)^ CLj- 



cl- 



■J-p- 



(18) 



However Eqs. (16,18) can only be consistent when both p and j have even values. Moreover 
the coefficients are independent of the sign of p and j. 
Then the pseudo-potential may be written as 

U m = E CL,p [D\^i + <2J [<;(/37) + D^iljih)] /(I + <^o.-)(l + V)- (19) 

-f'.lpl.b'l {even) 

according to the restrictions above. If limiting the summation to those terms with L equal 
to 2, one has 



U m = C2oo<o + 2C220<2 cos(27) I d^^(/?) 



(2) 



7(2)/ 



+2 



C2oo45 + 2C22242 cos(27)] dj^l/^) cos(27). 



(20) 



where (im^n(/3) is a reduced Wigner rotation matrix [14]. The second rank order parameters 
Oq q and (ig 2Cos(27) in this expression are related to the principal elements of the ordering 



matrix by 



Szz — d. 



(2) 
0,0 



(21) 



and 



SxX - Syy^VQ cos (27). 



(22) 



The order parameter rfg 2 cos(27) indicates the deviation of the ordering matrix from cylin- 
drical symmetry. For simplicity, it is convenient to write the pseudo-potential as 



C/(/37) = -kT a + b d'^m cos(27) I , 



(2), 



(23) 



where 



a — — 



C200 d'^ol + 2 C220 c?^2 cos (27) / /cT 



<2), 



(24) 



and 



6 = - 



2 C220 d^ol + 4 C222 <^ cos(27) /kT. 



(2), 



(25) 



The partition function for this single particle potential is then given by 

+ b dfl{(3) cos(27)] sin (3d(3 d-f. (26) 
The orientational molar potential energy and the molar entropy are given as follows 



.2. 

^ — / exp 
Jo Jo 



U = -{RT/2) adfl + h d^^l cos(27) 

a45 + ^43cos(27)| +i?ln(Z). 



5 = -i? 

respectively. Then the orientational molar free energy is straightforward: 

RT 



(27) 
(28) 



'%l + h 4,2 cos(27) -RT\n{Z). 

The order parameters are calculated from the equations below: 

/o " r <o(/?) exp [adgJ(/3) + hd^^l{(3) cos(27)] sin ^7 



(29) 



"0,0 — 



(30) 



and 



!T lo ^?2(/3) cos(27) exp [a<i(/3) + 6dg(/3) cos(27)l sin d^ 

dgcos(27) = ^ ^ . 



(31) 



It is worthwhile to note that d'^^ in Eq.(30) is simply P2, {P2 is the second Legendre polyno- 
mial). The free energy Fm will be a minimum provided the consistency equations for these 
order parameters are satisfied. With the aid of Eqs. (24,25), one can write the following 
equations: 



a — — 



a 0?^^ + 2A 4,2 cos(27) , 



(32) 



and 



b = 2Xa, 



(33) 



where 



a = -C2oo/kT. 



(34) 



A in Eq.(34) is the ratio C220/C200 a-nd it is a measure of the deviation from cyhndrical 
symmetry. If A is taken to be zero, the theory reduced to the Maier-Saupe mean field 
theory. 

Now we would hke to give the description of TsaUis thermostatistics with its axioms 
below. 



TsaUis thermostatistics (TT) has been extensively used to investigate the concepts of 
statistical mechanics and thermodynamics. In this context, it has been applied to various 
phenomena [17] and TT seems to be appropriate for the study of nonextensive systems. 
As mentioned above, we have choosen the nematic-isotropic transition as an application 
area of TT in earlier studies [11-13] and the present study would be the extensional study 
to investigate the nematic liquid crystals. In all of these earlier studies, we consider that 
the uniaxial nematic liquid crystals are formed by cylindrically symmetric molecules. How- 
ever, we wonder in this study that if the mesophase could be formed by non-cylindrically 
symmetric molecules, how does the nonextensivity affect the nematic-isotropic transition? 
Therefore our starting point is Luckhurst et al.'s theory [10] which is a molecular field the- 
ory and assumes the molecules of the nematic liquid crystals to be non-cylindrical and we 
examine the effects of the nonextensivity on the nematic-isotropic transition by considering 
this molecular field theory. 



B. Generalization of the theory 
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The understanding of the basic principles and properties of TT has gain fundamental 
importance. In this manner, TsaUis et al. have studied the role of the constraint within TT. 
Tsallis proposed an entropy definition in 1988 [18]: 

S,^k (35) 

where /c is a constant, Pi is the probabilty of the system in the i microstate, W is the total 
number of configurations and q is called the entropic index whose meaning will be given 
below. In the limit q ^ 1, the entropy defined in Eq.(35) reduced to the Boltzmann-Gibbs 
(BG) entropy. In other words, TT contains Boltzmann-Gibbs statistics as a special case 
with q = I. As well known, BG statistics is a powerful one to study a variety of the physical 
systems. However, it fails for the systems which (i) have long range interactions, {ii) have 
long range memory effects and (iii) evolve in multi-fractal space-time. The systems which 
has these properties is called "nonextensivc" and if a system does obey these restrictions, 
BG statistics seems to be inappropriate and one might need a nonextensive formalism to 
study on the physical system. 

It will be useful to write the pseudo-additivity entropy rule, 

S,(A + B)/k = (S,{A)/k) + {S,{B)/k) + (1 - ?) {S,{A)/k) {S,{B)/k) , (36) 

which reflects the character of the entropic index q and also of the nonextensivity. In this 
equation, A and B are two independent systems. Since all cases ^ (which is called 
nonnegativity property), g < 1, g = 1 and ? > 1, correspond to superadditivity , additivity 
and subadditivity respectively. 

Because of some unfamiliar consequences of first two energy constraints discussed in 
[13,19], Tsallis et al. present the third internal energy constraint as follows: 

^^^Mi = [/(3). (37) 

Z^i=l Pi 



This choice is commonly considered to study physical systems and is denoted as the Tsallis- 
Mendes-Plastino (TMP) choice. The optimization of Tsallis entropy given by Eq.(35) ac- 
cording to the third choice of the energy constraint results in 
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with 

1 

1-9 

(39) 

Therefore the ^'-expectation value of any observable is defined as 

{A^)g = (40) 

where A represents any observable quantity which commutes with Hamiltonian. As can be 
expected, when q — 1, the gr-expectation value of the observable reduces to the conventional 
one. At this stage, the important point is to solve Eq.(40) which is an implicit one, and 
Tsallis et al. suggest two different approaches; "iterative procedure" and —>■ /5" trans- 
formation. The nonextensivity, in fact, appears in systems where long range interactions 
and/or fractality exist, and such properties have been invoked in recent models of mangan- 
ites [20,21] as well as in the interpretation of experimental results. These properties are also 
appear in [22] where the role of the competition between different phases to the physical 
properties of these materials is emphasized. The formation of micro-clusters of competing 
phases, with fractal shapes, randomly distributed in the material is considered in [23,24], 
and the role of long range interactions to phase segregation [25,26]. 

If we look at the recent studies in which TT is employed, we see also that the entropic 
index is frequently established from the dynamics of the system under consideration an- 
alytically. Some of such studies are: Boghosian et al.'s study [27] for the incompressible 
Navier-Stokes equation in Boltzmann lattice models, Baldovin et al.'s studies [28,29] at 
the edge of chaos for the logistic map universality class, Oliveira et al.' studies [20,21] on 
manganites. 

According to these axioms of TT given above, the order parameters are calculated by 
using Eq.(40) instead of Eqs.(30) and (31). The exponential form of any function, say x, is 
given by 



w 



w 



4" = E 



exp (x) = [1 + (1 - q)x] 



(41) 



in Tsallis thermostatistics. Substituting this equation in the calculation, Eq.(40) is employed 
and then the self-consistent equations are solved to obtain the order parameters. 
The generalized free energy then follows within TT: 



(F) =^ 

{■L mjq 2 



a(d^ +b ((igcos(27))^ 



-RTlng{Zg), 



(42) 



where lnq{x) = {Z^ — 1)/(1 — q). 

Now it is the time to give the results of apphcation of TT to the Luckhurst et al.'s theory 
presented above. 



II. RESULTS AND DISCUSSION 

The first stage in our calculation is the identification of the transition from the uniax- 
ial nematic phase to the isotropic phase. This may be accomplished for a given molecular 
symmetry by determining the value of a at which the orientational free energy vanishes, 
provided there is no volume change at the transition. The results of such calculations are 
hsted Tables (1-3). Table 1 shows the transition temperatures, the critical values of order 
parameters (ig q and (ij) 2 cos(27) at the transition temperature and the entropy changes at the 
transition concerning A = 0.1 for some values of q. We can see immediately that increasing 
the entropic index increases ax and so decreases the nematic-isotropic transition temper- 
ature for constant C2oo- Similar behaviour is also seen from Tables (2) and (3) concerning 
A = 0.2 and 0.3 respectively. Also we see from Tables (1-3) that increasing the entropic 
index increases the order parameter ^\ at the transition for A = 0.1, 0.2 and 0.3. However 
there is a somewhat suprising result that as q is increased, the secondary order parameter 
(io^2cos(27) also increases at the transition for A = 0.2,0.3 while this is not the case for 
A = 0.1, that is, the secondary order parameter decreases with increasing qiox \— 1. We 
also observe from the Tables (1-3) that the entropy changes at the transition increases with 
increasing q. 
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In Fig.(l), we illustrate the dependence of the order parameter d^Q and on reduced 
temperature for some q values, with A = 0.1 only. The results are plotted as a function of the 
reduced variable {a/ax) which is identical to the reduced temperature T/T^ provided the 
coefficients C2jp are themselves independent of temperature. We observe that as q increased, 
the transition becomes more markedly first-order. It is well known from [10] that the effect 
of introducing deviations from cylindrical symmetry is to lower the curves for A = 0, 0.1, 0.2 
and 0.3 respectively, as well as changing their slopes. It is clear from Fig.(l) that as g is 
decreased, the similar behaviour is also seen, that is, the slope of the curve for A = 0.1 
decreases with increasing the entropic index q. The similar behaviour could also expected 
for other A values (A = 0, 0.1, 0.3). The secondary order parameter (iQ^2COs(27) is plotted as 

(2) 

a function of (Iqq, that is P2, in Fig. (2) with A = 0.2 for various values of the entropic index. 
These results exhibit an unusual behaviour, for the order parameter d!^^ cos (27) is observed 
to increase with increasing P2 pass through a maximum and then decrease. It is interesting 
to note that as q is decreased, the maximum value of the secondary order parameter assumes 
a higher value than those concerning lower value of q. 

Another important point is that Maier-Saupe theory gives a universal value of the ^q^q at 
the transition temperature {{P2)c = 0.429). However both Luckhurst et al.'s theory and the 
generalized form of the Maier-Saupe theory, including the present study, assume different 
values of this order parameter with changing A and q parameters. This fact is consistent 
with experiment in that experimental values of q for various nematics assume in a range of 
0.25 — 0.5 [30]. While A denotes the ratio C220/C20O) i-©- it is responsible from the deviations 
from the cylindrical symmetry, g is a measure of the nonextensivity of the system. But it 

is interesting that they exhibit some similar behaviour about the critical value of the order 

(2) 

parameter d^ q at the transition. 
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A = 0.1 

Q 




^ 2 




S iNk 


0.99 


4.458 


0.396 


0.01719 


0.348 


0.995 


4.481 


0.402 


0.01716 


0.357 


1 


4.505 


0.408 


0.0171 


0.384 


1.005 


4.528 


0.413 


0.017 


0.398 


1.01 


4.551 


0.417 


0.0169 


0.408 



Table 1. 



A = 0.2 
<l 


aK 


P2 


^Q,2 


Sc/Nk 


0.99 


4.335 


0.332 


0.0334 


0.256 


0.995 


4.36 


0.336 


0.0335 


0.265 


1 


4.386 


0.341 


0.0336 


0.275 


1.005 


4.412 


0.348 


0.0337 


0.301 


1.01 


4.438 


0.354 


0.0338 


0.314 



Table 2. 



A = 0.3 
<1 




P2 


^0,2 


SjNk 


0.99 


4.104 


0.2 


0.0399 


0.101 


0.995 


4.132 


0.203 


0.0404 


0.105 


1 


4.160 


0.207 


0.041 


0.112 


1.005 


4.1894 


0.212 


0.0417 


0.117 


1.01 


4.218 


0.218 


0.0424 


0.126 



Table 3. 



III. CONCLUSION 

Many molecular theories assume that an essential characteristic of compounds forming 
liquid crystals is the rod-like shape of their constituent molecules that means an high length 
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to breadth ratio. This assumption leads to the approximation that the molecules might be 
assumed to have cylindrical symmetry. However the molecules, in fact, are lath-like and thus 
do not possess the high symmetry. In this sense, a molecular theory was developed [10] for 
an ensemble of such particles based upon a general expansion of the pairwise intermolecular 
potential together with the molecular field approximation. 

Tsallis thermostatistics has been commonly employed in studying on the physical systems 
as a nonextensive statistics and many studies are investigating the nonextensive effects 
in different systems and phenomena. Thus we would like to handle, in this study, this 
molecular field theory [10] and to investigate the effects of the nonextensivity for uniaxial 
nematics formed by non-cylindrically symmetric molecules. With this aim, we investigate 
the dependence of the long range order parameter on reduced temperature and report the 
variation of the critical values of the order parameters and entropy change at the transition 
temperature with the entropic index. 
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IV. FIGURE AND TABLE CAPTIONS 



Table 1. The order parameters and entropy change at the nematic-isotropic transition 
with A = 0.1 for some q values. 

Table 2. The order parameters and entropy change at the nematic-isotropic transition 
with A = 0.2 for some q values. 

Table 3. The order parameters and entropy change at the nematic-isotropic transition 
with A = 0.3 for some q values. 

Figure 1. The variation of the order parameters (Tq q and ^0,2^03(27) with the reduced 
temperature for some q values, with A = 0.2. q = 1 represents the Luckhurst et al.'s theory. 
Figure 2. The dependence of do^2cos(27) on d!ifi{— P2) for various q values. 
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